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The study of non- Abelian Majorana zero modes 
advances our understanding of the fundamental 
physics in quantum matter, and pushes their po- 
tential applications to topological quantum com- 
putation. Recent investigations have shown that 
in the two-dimensional (2D) and ID chiral su- 
perconductors, braiding isolated Majorana zero 
modes leads to noncommutative transformation 
of the quantum states [UlS]. However, Majorana 
fermions in a Z2 time-reversal invariant (TRI) 
topological superconductor come in pairs due to 
Kramers' theorem [3]. Therefore, braiding oper- 
ations in TRI superconductors always exchange 
two pairs of Majorana fermions. In this work, 
we show surprisingly that, due to the protection 
of time-reversal symmetry, non- Abelian statistics 
can be obtained in the ID TRI topological su- 
perconductor and may have advantages in apply- 
ing to topological quantum computation. Fur- 
thermore, we predict a novel phenomenon in the 
Josephson effect for ID TRI topological super- 
conductors, which measures the topological qubit 
states in such systems. 

The search for exotic non-AbeUan quasiparticles has 
been a focus of both theoretical and experimental stud- 
ies in condensed matter physics, driven largely by their 
promising applications to fault-tolerant topological quan- 
tum computation [3H3- Following this pursuit, the topo- 
logical superconductors have been brought to the fore- 
front for they host exotic zero energy states known as 
Majorana fermions p!U1115j . For 2D chiral p + ip pairing 
state, which breaks time-reversal symmetry, one Majo- 
rana mode exists in each vortex core ^, and for ID p- 
wave case, such state locates at each end of the system 
[5]. Existence of 2n Majorana zero modes leads to 2"~^- 
fold ground-state degeneracy, and braiding two of such 
modes transforms one state into another which defines 
the non- Abelian statistics [HIS]- Remarkably, Majorana 
end states have been suggestively observed through tun- 
neling measurements |16j in ID effective p-wave super- 
conductors obtained by semiconductor nanowire/s-wave 
superconductor heterostructures [rTHTO] 

Recently, a new class of topological superconductors 
with time-reversal symmetry, referred to as Dili symme- 
try class superconductor and classified by Z2 topological 
invariant [21 [2011231 . have attracted rapidly growing ef- 
forts [24ll28) . Different from chiral superconductors, in 
Dili class superconductor the zero modes come in pairs 
due to Krammers' theorem. Several interesting proposals 



have been studied to realize Z2 TRI Majorana quantum 
wires using proximity effects of d-wave, p-wave or si- 
wave superconductors. It was shown that at each end of 
such a quantum wire are localized two Majorana fermions 
which form a Kramers doublet and are protected by time- 
reversal symmetry [2S1I28) . 

With the practicability in realization, a fundamental 
question is that can the Dili class topological super- 
conductor be applied to topological quantum computa- 
tion? The puzzle arises from the fact that braiding the 
end states in a Dili class ID superconductor always ex- 
changes two pairs of Majorana modes as illustrated in 
Fig. [1] While braiding two pairs of Majoranas in chiral 
topological superconductors yields Abelian operations, in 
this work, we show that braiding Majorana end states in 
Dili class superconductors is non- Abelian due to the pro- 
tection by time-reversal symmetry. 
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FIG. 1: Braiding Majorana zero modes in ID TRI 
topological superconductor through a T-junction fol- 
lowing the study by Alicea and the colleagues [2j. 
a, The Dili class topological superconductor can be realized 
by depositing a conducting nanowire (NW) on the top of a 
non-centrosymmetric superconductor (SC). b-e. Braiding end 
modes through a T-junction. The dark (light gray) area of 
the nanowires depicts the topological (trivial) region, which is 
controllable by tuning the chemical potential in the nanowire. 
The arrows depict the direction that the Majorana fermions 
are transported to in the braiding process. 

In the Methods section, we show that a ID TRI topo- 
logical superconductor can be engineered by inducing p- 
wave superconductivity in a conducting wire in proximity 
to a non-centrosymmetric superconductor (Fig. flk) 29] . 
In the topological phase, at each end of the wire are local- 
ized two Majorana zero modes 7^ and 7^ (j — 1,2), trans- 
formed by time-reversal operator that T^^^jT — 7j and 
T^^"ijT = — 7j [20] . Note that braiding Majorana end 
modes is not well-defined for a single ID nanowire and, as 



first recognized by Alicea et al., the minimum setup for 
braiding requires a trijunction, e.g. a T-junction com- 
posed of two nanowire segments [5]. The braiding can 
be performed by transporting the Majorana zero modes 
foUowing the steps as iUustrated in Fig. [ijb-e) . 

To visuahze the braiding rule, let us consider first the 
simplest situation that the Dili class topological super- 
conductor is composed of two decoupled copies (e.g. cor- 
responding to spin-up and spin-down, respectively) of 
ID chiral p-wave superconductors. It is then straight- 
forward to know that the whole braiding is a product 
of two independent processes of braiding 71^2 and 71^2, 
respectively, which gives the TRI unitary transforma- 
tion matrix Ui2{T,f) = exp(|7i72) exp(|7i72) follow- 
ing the studies by Ivanov, Alicea and the colleagues 
^ 12]. This braiding is nontrivial as presented below, 
and detailed in the Supplementary Material. We con- 
sider two Dili class wires with eight Majorana modes 
7i,...,4 and 71,..., 4 (Fig. uh), which define four complex 
fermion modes by /i = 5(71 + il2),f2 = \{lz + ili), 
and /i^2 — T~^fi.2T- The Hilbert space of the four 
complex fermions is spanned by sixteen qubit states 
\nin2)\n,h2) = (/J)"H/1)"^ (/7)"^ (/|)"=|00)|00). If the 
initial state of the system is |00)L|00)ii:, for instance, by 
braiding the two pairs of Majoranas 72, 72 and 73, 73 we 
get straightforwardly 

C/23(T,f)|00)L|00)fl, = ^{\06)LmR+\li)L\ll)R 

+ i|io)L|io)fl-z|oi)i|oi);^)(i) 

Here L/R represents the left/right nanowire segment. It 
is interesting that the above state is generically a four- 
particle entangled state, which shows the natural advan- 
tage in generating multi-particle entangled state using 
Dili class topological superconductors. Furthermore, a 
full braiding, i.e. braiding twice 72,72 and 73,73 yields 
the final state |11)l|11)_r, which distinguishes from the 
initial state in that each copy of the p-wave superconduc- 
tor changes fermion parity. 

From the above discussion we find that in the braiding 
the Majorana modes 7^ do not feel their time-reversal 
partners jj , which is an essential difference from the sit- 
uation in exchanging two pairs of Majorana fermions in a 
chiral superconductor, and makes the braiding operator 
in a TRI topological superconductor nontrivial (Fig. ^p- 
c). As presented below, the above braiding statistics, 
while being obtained in the special situation, are valid 
in the generic Dili class ID topological superconductors 
and protected by time-reversal symmetry. 

We first prove a central result that by grouping all the 
quasiparticle states into two sectors being time-reversed 
partners of each other, the fermion parity is conserved 
for each sector, not only for the total system. The 
proof is equivalent to showing that in a TRI Majo- 
rana quantum wire, the four topological qubit states 
jnifii) (ni, fii = 0, 1) are decoupled from each other with 
the presence of finite TRI perturbations. The coupling 
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FIG. 2: Non-Abelian statistics in Dili class ID topo- 
logical superconductor, a, Majorana end modes 72, 72 and 
73,73 are braided tlirough similar processes shown in Fig. [T] 
b, Braiding Majorana modes in Dili class superconductor is 
equivalent to two independent processes of exchanging 72,73 
and 72, 73, respectively. In the depicted process 72 (72) crosses 
only the branch cut of 73 (73), and therefore acquires a minus 
sign after braiding, c, In contrast, if braiding two Majorana 
pairs in a chiral superconductor, for the depicted process 72 
(72) crosses the branch cuts of both 73 and 73, and then 
no sign change occurs after braiding 1^. Therefore, braiding 
twice two Majorana pairs always returns to the original state. 



Hamiltonian, assumed to depend on a manipulatable pa- 
rameter A, should take the generic TRI form V{X) = 
i£^i(A)(7i72 - 7172) + i£'2(A)(7i72 - 727i)> which splits 
the two even parity eigenstates |00) and |11) by an energy 
E{X) 



2^yEf + El Since |10) and |0i) form a Kramers' 
doublet, the transition between them is forbidden by 
time-reversal symmetry. Then the fermion parity con- 
servation requires that the following adiabatic condition 
be satisfied in the manipulation: |(11|A(9a|00)| ^ 2E{\), 
where A = dX/dt. This is followed by 



R: 



,dXde, 



2E{X) I dt dX 



< 1, 6* = tan" 



E\ 
E2 



(2) 



In a realistic system, the coupling £^1^2 are controlled 
by adjusting the bulk gap i^^buik or the distance d be- 
tween the Majorana modes, and are generically exponen- 
tial functions of E'buik'^- Since 7^ and 7^ are connected 
by T-transformation, their wave functions have exactly 
the same spatial profile, which leads to the same expo- 
nential form of the couplings i?i,2(A) = ai^2{X)e~^^^'^^^'^ . 
The pre-factors aj[X) depend on the local couplings (e.g. 
hopping terms, pairings, etc.) between electrons belong- 
ing to the same (for j — 1) or different (for j = 2) sectors 
of the time-reversal partners. For constant and homoge- 
neous local couplings, we expect that while the magni- 
tudes of otj{X) depend on A, their ratio axja^ is nearly a 
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FIG. 3: Adiabatic condition and fermion parity conser- 
vation for each sector of time-reversal partners, a-b, 

The couplings £1,2 between Majorana end modes are manip- 
ulated by tuning the chemical which changes the bulk gap 
(A = i?buik) in the nanowire (a), and by varying the length 
(A = d) of a trivial region (gray color) which separates the 
two pairs of Majorana modes (b). c-d, The energy splitting 
E between |00) and |11) (red curves) and the ratio R (blue 
curves), as functions of iSbuik (c), and versus the trivial re- 
gion distance d (d). The parameters in the nanowire are taken 
that the proximity induced p-wave pairing Ap — l.OmeV, s- 
wave pairing Aa — 0.5meV, and the spin-orbit coupling en- 
ergy Eso ~ O.lmeV. In the numerical simulation we assume 
that the coupling energy E is tuned from to l.OmeV in 
the time l.Ofis. We also numerically confirmed the adiabatic 
condition R <ti 1 with other different parameter regimes. 



constant (rigorous proof is given in Supplementary Mate- 
rial). Therefore we always have d\9 ~ 0, which validates 
the adiabatic condition. This general statement is clearly 
shown with numerical results in Fig. [3] The fermion par- 
ity conservation for each sector shows that an isolated 
Dili class ID Majorana wire should stay in one of the 
four fermion parity eigenstates germinated by non-local 
complex fermion operators fj and fj, given that time- 
reversal symmetry is not broken. In particular, one can 
always prepare a nanowire initially in the ground state 
jOO) or 1 11) by controlling the initial couplings £'1.2 (A), 
and then manipulate the states adiabatically. 

The fermion parity conservation for each sector im- 
plies that the exchange of Majorana end modes in Dili 
class superconductor is generically equivalent to two in- 
dependent processes of braiding Majorana fermions of 
two different sectors, respectively. This is because, first 
of all, braiding adiabatically the Majorana zero modes, 
e.g. 71, 71 and 72,72, does not affect the bulk states 
which are gapped. Furthermore, assuming that other 
Majorana zero modes are located far away from 71,2 and 
7i 2, the braiding evolves only the Majorana modes which 
are exchanged. Finally, due to the fermion parity conser- 



vation, in the braiding the fermion modes /i and /i are 
always decoupled and their dynamics can be derived in- 
dependently. We therefore obtain the braiding matrix by 
Ui2{T,f) = exp(|7i72) exp(|7i72), as presented above. 
A rigorous derivative of the braiding matrix is given in 
the Supplementary Material. 

It is important to study how to detect the topological 
qubit states in a Dili class Majorana quantum wire. We 
consider a Josephson junction shown in Fig. Hk formed 
by Dili class superconductor. As derived in the Supple- 
mentary Material, the effective coupling Hamiltonian of 
the Josephson junction is given by 



^ Josephson junction " 

(l)=0 ]r ^^=^ ^ 





FIG. 4: Josephson measurement of the topological 
qubit states in Dili class ID topological supercon- 
ductor, a, The sketch of a Josephson junction with phase 
difference <j). b, The single particle Andreev bound state spec- 
tra versus the phase difference cj>. c, The energy spectra of the 
four qubit states |nini) (ni, h\ — 0, 1) according to the results 
in b. d, The Josephson currents (in units of 2eAp/h) for dif- 
ferent topological qubit states. Parameters used in the numer- 
ical calculation are taken that Ap = l.OmeV, As = 0.25meV, 
Eso ~ O.lmeV, the width of the junction d — 0.5^, and in 
middle trivial region (gray color) of the junction the chemical 
potential is set to be at the band bottom. 



^off — 



iTo cos -(7172 - 7172) + i^i sin0(7i7i - 7272)(3) 



where (j) is the phase difference across the junction. The 
Fo-tcrm in Heff represents the first-order direct coupling 
between Majorana fermions at different junction leads. 
The Ti-term is resulted from the second-order perturba- 
tion of the tunneling process, and this term vanishes if 
the s-wave pairing A^ = (a detailed discussion can be 
found in the Supplementary Material). Redefining the 
Majorana bases by 7i = 71 + 72 and 72 = 72+71, we 
recast the above Hamiltonian into Hcs — J(roCOS(/)/2 -|- 
Fi sin 0)7^ 72 — i(FoCOS(/)/2 — Fi sin</))7^72, which gives 
straightforwardly the Andreev bound state spectra. Nu- 
merical results are shown in Fig. |41d-c. The Josephson 



currents are obtained by the slope of the Andreev bound 
state spectra. In particular, we have that the Josephson 
currents JT""^ = ±|rosin^ for the even parity states, 

|06) and 1 11), and Jf'^ = ±f Fi cos0 for the odd parity 
states, |0i) and |10), respectively (Fig. |4|i). 

Interestingly, the currents for odd-parity states are 
of 27r periodicity, half of those for even-parity states 
(Fig. Eli). This reflects that J|™" is contributed from 
the dnect Majorana coupling induced by first-order tun- 
neling, while JS'^'^ is a consequence of the second-order 
tunneling process which corresponds to the Cooper pair 
tunneling. Such result is also consistent with the fact that 
the time-reversal symmetry is restored with |01) and |10) 
forming Kramers' doublet at cj) — mn, which necessitates 
the 27r periodicity in their spectra. Furthermore, the two 
qubit states with the same total parity (e.g. |00) and 11)) 
are distinguished by the direction of the currents. The 
qualitative difference in the Josephson currents provides 
direct measurements of the four topological qubit states 
in the experiment. 



Methods 

Topological superconductor of Dili class by 
proximity effect. Several interesting proposals have 
been studied for realizing Dili class ID topological su- 
perconductor |26H28j . Here we consider to engineer such 
topological superconductor by depositing a conducting 
quantum wire on a non-centrosymmetric superconduc- 
tor (Fig. II]) . Details of the proximity effect are given in 
the Supplementary Material. The induced s- and p-wave 
superconductivity in the wire can be obtained by inte- 
grating out the degree of freedom of the superconductor 
substrate. This gives the spectral function of the wire 



A{UJ,k,j;) 



2tt 



r{tr[T2G'wirc(w + «0+, /C^r)]} 



where G^^.^^uj, k^) = ujI-^k^T^ - S(a;, k^) is the Green's 
function of the wire. The induced s- and p-wave pairings. 
As and Ap, are encoded in the self-energy Yi{u],kx) — 
|i_Lp / dky[(jj — Ha,c{krc, ky)]~-^ duc to the coupling to the 
bulk superconductor with its Hamiltonian denoted by 
Hsc- Here t± characterizes the coupling between the 
nanowire and substrate superconductor through their in- 
terface. The spectral function determines the effective 
Hamiltonian and the band structure of the wire. In the 
Supplementary Material we find that that the nanowire 
is in the topological regime when its chemical potential 
/i the induced pairings satisfy —2t^ < fi < 2t^ and 



|Ap| > I As I, where t^ denotes the hopping coefficient 
in the nanowire. 

Non-Abelian braiding. To realize a Dili class su- 
perconductor applies no external magnetic field, which 
can be advantageous to construct realistic Majorana net- 
work to implement braiding operations. In compari- 
son, for the chiral topological superconductor observed 
in a spin-orbit coupled semiconductor nanowire using s- 
wave superconducting proximity effect pTHTO) . the ex- 
ternal magnetic field should be applied perpendicular 
to the spin quantization axis by spin-orbit interaction, 
driving optimally the nanowire into topological phase 
[T71 [ini [30] ■ It is shown that for a network formed by mul- 
tiple nanowire segments, such optimal condition cannot 
be reached for all segments without inducing detrimental 
orbital effects, which creates further experimental chal- 
lenges in braiding Majoranas ISO]. It is clear that such 
intrinsic difficulty is absent in the present Dili class TRI 
topological superconductor, and one may have more flexi- 
bility in constructing 2D and even 3D Majorana networks 
for topological quantum computation. 

Josephson effect in Dili class topological su- 
perconductor. After braiding the end modes, we can 
transport the Majorana modes to the Josephson junction 
and measure the topological qubit states by Josephson 
effect. It turns out that the first-order direct coupling 
occurs only between Majorana modes at different ends 
of the junction, yielding the first term in the coupling 
Hamiltonian Hcs of Eq. (|3|. The leading-order contri- 
bution to the coupling between Majorana fermions at 
the same end (ijjjj) comes from the second-order per- 
turbation in the tunneling process. This is because, the 
coupling term such as ijj'jj breaks time-reversal symme- 
try, while the direct coupling between jj and 7^ does not 
experience the phase difference across the junction and 
therefore should preserve time-reversal symmetry. To get 
the second term in Hcs, the minimum requirement is 
to consider the second-order tunneling process. Further- 
more, since the system restores time-reversal symmetry 
at = TTiTT, the leading-order contribution to i"fj"fj must 
be proportional to sin 0, which has 2tt periodicity. These 
properties are unaffected by TRI disorder scattering. 
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Supplementary Information 

In this Supplementary Material we provide the details of the results presented in the main text. In particular, we 
study in detail the proximity effect of a non-centrosymmetric superconductor in a conducting nanowire, the fermion 
parity conservation for each sector of the time-reversal partners, the braiding statistics of Majorana end modes in 
Dili class ID topological superconductor, and the Josephson effect for the measurement of topological qubit states in 
a Dili class superconductor. 

S-1. TOPOLOGICAL SUPERCONDUCTOR OF Dili CLASS BY PROXIMITY EFFECT 



The Dili class ID topological superconductor can be engineered by depositing a conducting quantum wire on a non- 
centrosymmetric superconductor thin-film which induces s- and p-wave pairings in the wire by proximity effect. The 
total Hamiltonian of the heterostructure system reads H = Hgc + H.^i^g + H^, where 7?sC:^wirc: and H^ represent the 
Hamiltonian for the substrate superconductor, the conducting wire, and the tunneling at the interface, respectively. 



Due to the presence of spin-orbit coupling, a typical non-centrosymmetric superconductor has both p-wave and s-wave 
pairings. To differentiate from the notations used for the nanowire, we denote by the pairings in the superconductor 
by As ' and Ap ' , respectively. The BdG Hamiltonian for the 2D non-centrosymmetric superconductor is given by 



Hsc = ^ [£.{kx, ky)Tz + a^j^' sin kya^ - a)^ sin k^ayTz + A^°' sin kyCr^Tx - A^"^ sin k^Ty -~ Af>ayTy] , (SI) 

where S,{kx , ky) — ~2t'^^^ {cos k^ + cos ky) — ^j,^"^ is the normal dispersion relation with i^^^ the hopping coefhcient in the 
superconductor, aj and Tj [j — x, y, z) are the Pauli matrices acting on the spin and Nambu spaces, respectively, a}j 
is the spin-orbit coupling coefficient, and /i'"' is chemical potential. The pairing order parameters can be reorganized 
by A = {As -f- d • a){iay), with the d- vectors defined as d = Ap '(— sinfey,sinfc2;,0). In the above Hamiltonian we 
take that the g- vector for spin-orbit coupling, given by g = a)^ {sin ky, — sinfe^;), is antiparallel to d. 

A single-channel ID conducting quantum wire, being put along the x axis, can be described by the following 
Hamiltonian 






^wirc = > (-'2tu, COskx - fly,)T^, (S2) 






with t„, the hopping coefficient and fi^ the chemical potential in the wire. It is noteworthy that the intrinsic spin-orbit 
interaction is not needed to reach the TRI topological superconducting phase in the nanowire, while the proximity 
effect can induce an effective spin-orbit interaction in the wire. Finally, we give the tunneling Hamiltonian Ht for 
the interface. Assuming that at the interface the coupling between the substrate superconductor and the nanowire 
is uniform along the wire, the momentum k^ is still a good quantum number (we note that this assumption can be 
relaxed for the existence of Majorana fermions in the topological phase). Then the tunneling Hamiltonian can be 
written down as 

Ht = -ti_ Y. clJk,)d^,Ak.) + h.c, (S3) 

k^ ,(7 

where a ="[, J, are the spin indices, t± denotes the tunneling coefficient between the nanowire and substrate supercon- 
ductor, cl ^,Ci ,a and d\ ^,di ^^ are the creation and annihilation operators of electrons for the quantum wire and 
the superconductor, respectively. The site number iy characterizes where the heterostructure locates along y axis in 
the non-centrosymmetric superconductor. 

The induced superconductivity in the wire can be obtained by integrating out the degree of freedom of the su- 
perconductor substrate, and will be calculated numerically, as described below. We perform the integration in two 
steps. First, consider that the non-centrosymmetric superconductor is uniform, and we determine the Green's function 
Gs{kx,iy) of the substrate superconductor with momentum kx and at the site iy below the nanowire by standard 
recursive method [1]. Then, the coupling of the nanowire to the superconductor can be reduced to the coupling to 
the site iy below the wire and described by the Green's function Gs{kx,iy)- Integrating out the degree of freedom of 
the sites in the superconductor right below the nanowire yields a self-energy for the Green's function of the nanowire, 
which gives rise to the proximity effect. The effective Green's function of the nanowire takes the form 

Gll,{iu;, kx) = -. ^ ^,^-y (S4) 

where the self-energy reads S(iix)) = t\Gs{kx,iy)- Finally the spectral function is determined by 

A{i^.k,) = -^3{Tr[T,G;f,,(c. + *0+,fc,)]}, (S5) 

with 5 taking the imaginary part, Tr denoting the trace over the spin and Nambu spaces, and 0"*" a positive infinites- 
imal. The spectral function determines the bulk band structure, which is numerically shown in Fig. [ST] with different 
chemical potentials of the nanowire. In particular, from the numerical results we find that the nanowire is in the 
topologically nontrivial regime when |/iu,| < 2|i^| and |Ap '| > |Ai '| which leads to the induced pairings in the wire 
|Ap| > |As|, while it is in the trivial regime when |Ai | > |Ap | or |/i„| > 2|i^| (i.e. the chemical potential is tuned 
out of the band of the wire). When tuning the chemical potential down to the band bottom, the bulk gap in the 
nanowire is reduced and closes right at the bottom, implying the critical value of the chemical potential /^^ = — 2i^ 



(similar results can be obtained around 2tu,, the top of the band) [Fig. SI (a-b)]. In the topological regime at each 
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FIG. SI: The logarithmic plot of the spectral function for the nanowire/non-centrosymmetric superconductor heterostructure. 
The yellow dotted curves show the bulk band structure of the nanowire system. The red solid areas (in the upper and lower 
positions of each panel) represent the bulk states of the substrate superconductor, (a) Topological regime with the chemical 
pote ntia l /i- 
[Fig 



S2 
point fj,^ 



set as — 2t„ + 5|Ap | in the nanowire. In this regime at each end of the wire localized two Majorana zero modes 



(b) Topological regime with reduced bulk gap by tuning /i^ 



2tw + |Ap I close to the band bottom, (c) Critical 



— —2tn] for the topological phase transition with the bulk gap closed, (d) Trivial phase regime for the nanowire with 



/i — —2tm — |Ap |. Other parameters are taken that t± — 0.5tiu — O.St'"', jAs| = 0.5|Ap|, and q)j = |Ap| 



end of the nanowire are localized two Majorana zero modes "fj and jj (j — L, R), with their wave functions shown in 
Fig. |S2| Further lowering the chemical potential reopens the bulk gap, and the system is driven into a trivial phase 



[Fig. SI (d)] 



It is interesting that the phase diagram in the nanowire does not depend on the parameter details of the couplings 
between the wire and the substrate superconductor, and the topological regime in the wire can be obtained in a large 
parameter range that —2tw < Hw < '2tw when jAJ, | > jAsI*-"-'. With such a large window for the topological phase, 
the present proposal shows advantages in engineering the Dili class topological states of the wire by simply tuning 
fiw to be below or above the band bottom of the wire. 



S-2. FERMI PARITY CONSERVATION 



Fermion parity measures the even and odd numbers of the fermions in a quantum system. For a superconductor 
because of the pairing gap the number of fermions can only vary by pairs, which leads to the fermion parity conserva- 
tion. Therefore, the total fermion parity is always conserved in a Dili class ID topological superconductor. Moreover, 
in this section we show in detail that in a Dili class ID topological SC the fermion parity is actually conserved for 
each sector of the time-reversal partners, more than for the total system. It is trivial to know that this result is true if 
the Dili class topological superconductor is composed of two decoupled copies (for instance, corresponding to spin-up 
and spin-down, respectively) of ID chiral p-wave superconductors. The each copy of the chiral p-wave superconductor 
conserves the fermion parity conservation. Therefore, for the proof we shall focus on the generic Dili class topological 
superconductor where the TRI couplings are allowed between different copies of time-reversal partners. 



(a) 




(b) 0.5 
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FIG. S2: (a) Two Majorana bound modes exist at each end of the nanowire in the topological regime with fj, = —2tm + \Ap \ 
(b-c) The wave functions of two Majorana modes "/l/r and Jl/r at the same end have exactly the 
1 ^ the coherence length in the nanowire. 



as considered in Fig. SI 



same spatial profile, wit 



We consider a single Majorana quantum wire, which hosts four Majorana end modes denoted by 71 2 and 71^2 £^nd 
transformed by time-reversal operator that T~^7jT = jj and T^^jjT — —jj |2J. With the four Majorana states we 
can define two non-local complex fermions by /i — |(7i -|- 1^2), fi = 5(71 — «72), which germinate four topological 
qubit states \nihi) with ni,ni = 0,1. The proof of fermion parity conservation for each sector is equivalent to 
showing that the four topological qubit states \nihi) are generically decoupled from each other with the presence of 
TRI perturbations in the Majorana modes. 

Note that the coupling between the Majorana modes localized at the same end of the nanowire, jj and jj, breaks 
time-reversal symmetry. The coupling Hamiltonian in terms of Majorana end modes should take the following generic 
TRI form 



V{X) = i^i(A)(7i72 - 7172) + i£^2(A)(7i72 - 7271), 



(S6) 



where we assume that the couplings £'1,2 (A) depend on an experimentally manipulatable parameter A. The above 
Hamiltonian can be rewritten in the block diagonal form with new Majorana bases that 



y(A) = *ii;(A) [7(^)7(2) -7(1)7(2)] 



(S7) 



^(1) = ^, -^(2) ^ 



where 7(1) = 

is defined via ta.n6 = Ei/E2. The complex fermions /(i) and /(i) in the eigen-basis are then defined by 



71 7 7^"' = 7117^"' — sin 6*72 -f 008^72, 7(2) — sin072 — cos 6*72 , and E — y/Ef + _E| . The mixing angle 6 



f'' = '^H'' 



Z7 



(2)1 



/(I) = i[7(i)-z7(2)] 



(S8) 
(S9) 



It is easy to know that the even parity eigenstates |00) and |li) germinated by /(i) and /(i) acquire an energy splitting 
2_E(A), while the odd parity states |01) and |10) are still degenerate. To prove the fermion parity conservation for each 
sector, we need to confirm that all the four topological qubit states \nini) can evolve adiabatically when the coupling 
Hamiltonian V{\) changes with the parameter A. Since |10) and |01) form a Kramers' doublet, the transition between 
them is forbidden by the time-reversal symmetry. Therefore, we only need to consider the adiabatic condition for the 
two even parity states. The fermion parity conservation for each sector is guaranteed when the following adiabatic 
condition is satisfied in the manipulation 



(ll|--|00)|«2|i.(A)|. 



(SIO) 



9 

It should be noted that the adiabatic condition needs to be justified only in the presence of finite couplings. When 
E{X) — >■ 0, the couplings between Majorana end modes vanish and then all the topological qubit states are automati- 
cally decoupled from each other. One can verify that 

In the above formulas the derivative of the bases jj , jj with respect to A does not contribute to the left hand side of 



Eq. (SIO), and is neglected. The condition (SIO) then reads 

|f|^|«4|^(A)|. (SI3) 

We show below that the above condition is generically satisfied in the realistic materials. 

According to the the previous section, with the proximity induced p-wave and s-wave superconducting pairings, the 
effective tight-binding Hamiltonian in the nanowire can be written as 

where the hopping coefhcients and the chemical potential are generically renormalized by the proximity effect. Without 
loss of generality, in the above Hamiltonian we have taken into account the spin-orbit interaction described by the 
tf- term, and the random on-site disorder potential V.'^'*' with (K'^'^) — 0. For the case with uniform pairing orders, 
the parameters A^ and Ap can be taken as real. On the other hand, for one-dimensional system, the phases in the 
(spin-orbit) hopping coefficients can always be absorbed into electron operators. Therefore, in the following study we 
consider that all the parameters in H^^.^ are real numbers. 

In the topological regime, at each end of the wire we obtain two Majorana zero modes which are transformed to 
each other by time-reversal operator. In terms of the electron operators, these bound modes take the form 

j 

3 

72 = i^[uf\xj)cfixj)+uf\xj)c^{xj)~-u'f^*{xj)c\{xj)~-uf^*{xj)cl{xj)], (S17) 

j 
72 = i^[uf''*{xj)c^{xj)-uf^*{xj)c-f{xj)-uf\x^)cl{xj)+uf\xj)c\{xj)]. (S18) 



Since the coefficients in H^^^ are real, we have that u[]' — u[]'* . The coupling energies between the Majorana 



dc 

*(7i|-f^wfrol72) = -«(7i|-f^wLl72)- Using the relations 



modes a,t left (71,71) and right (72,72) ends are calculated by Ei = i{^i\H^%Jj2) = -i(7i|i?°frcl72) and E2 



Co 



,(^h^\^.^ _„(i)^^.>.~. _,-„(2)/^.^__L,-„(2) 



■j^ ~ uy{xj)^i-uy{xj)^i-iuy{xj)"f2 + iul'{xj)^2, (S19) 



c. 



M)f^ ^„, ^,,(1)^^ ^^, ,-,,(2)/^ \^, ,-„(2) 



'jl - "1 i^jhi + ""I {^j)li - *4 (a;j)72 - iu\ '{xj)^2, (S20) 



10 
we obtain that 

+ Y,A,[u\'\x,)uf{x,)~u['\x,)uf\x,)]+Y,Vf'H'\x,W^Hx,), (S21) 

+ E t7^H'^i^^)u?\^J) ~ v^f'ix^v^ix^ + E ^s{uf{x,)uf{x,) + ^^^ (x,- )4^) (:.,)] 

+ E^''1"r(^^>r(^^)-"r(^^-)4'H^.-)- (S22) 

Note that tf° = — tj° and for a uniform nanowire we expect that ^,^ ■> ui ^(a;^)^^ (a;^) = ^/^ v wi '(xi)M| (a;^) and 
^ wi ^(a;j)w| (a;j) = ^ ■ ui '(xj)u| (xj). With these properties we find that in Ex the terms corresponding to if° 
and As vanish, while in E^ the terms for i^j, Ap, and V,!^"^ vanish. We then have 

i?i = E U,u^^\^^)u^^\x,)+ E Af^.^.«(x,)^.(2)(a;,.)+EV;'M'H^.)4''(^,), (S23) 

The wave functions of Majorana bound modes decay exponentiahy as a function of the distance from the end of 
the nanowire, muhiplying by an oscillatory function with the oscillating period equal to the Fermi wavelength in 
the nanowire. This implies that uh oc sin(fcp'a:)e^^/'» and uh oc sin[A;j;-(L — a;)]e^(^^^)/^, where ^ is the effective 
coherence length of the wire. In the realistic material, we consider that the chemical potential in the nanowire is far 
below the half-filling condition and thus kpa <^ 1. In this way we have uh {xj)ua (xj) ~ Ua {xj)ua {xj±i)e^°''^. 
Furthermore, the coherence length is typically much larger than the lattice constant ^ S> a, and we can further 
approximate that m^ {xj)ucr {xj) ~ Ua {xj)ua- {xj±i). Bearing this result in mind we get 

Ex ^ J2 uA'H^^W'\^,)+ E A^A'^ix,)u^^Hx,) + J2vf^^ui'H^j>^'H^j), (S25) 

E^ = Y. tT,u'^^\x,)u'^^\x,)+ E A,z.(i)(x,)z.(2)(a;,). (S26) 

The spin-orbit hopping coefficient t^° = — i^° and the p-wave pairing A^ = — A^^ are staggered parameters. In the 
limit that kpa <^ 1 and ^ 3> a, the summation for such two terms in Ei and E2 also turns out to be zero. On 
the other hand, the spin-conserved hopping is a constant and we denote t^- — tji — t. Finally, if the the random 
potential V,'^'^ with (V^,'*"') = is distributed homogeneously in the nanowire, we expect that the last term in Ei gives 

Vb X]i cr "o- {^j)'^<y i^j) with the constant factor Vq depending on the specific disorder profile and much less than the 
amplitude of the disorder potential. The couplings _Ei 2 become 

Ex ^ (t + Vo) E ui'H^^)u^aH^J), (S27) 

S2 ^ A, E U^a\^^W'\x,)■ (S28) 

From the above result we find that E2/EX ~ As/(i + Vq), which is consistent with the fact that when A^ = the 



original Hamiltonian (S14| can be block diagonalized and then E2 — 0. This implies that in the realistic nanowire 
materials while the magnitudes of i?i,2(A) depend on A which determines the overlapping between the wave functions 
of Majorana bound modes at left and right ends, their ratio is nearly a constant. Therefore we always have 

|.0, ,S29, 
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which vahdates the adiabatic condition. The results in the Eqs. (S27) and (S28l have a simple physical picture. 



Being proportional to the overlapping between the wave functions of Majorana bound modes at different ends, the 
couplings El 2 are exponential decaying functions of the nanowire length. Since the Majorana modes jj and 7j are 
connected by T-transformation, their wave functions have exactly the same spatial profile, which leads to the same 
exponential form for the couplings i?i.2(A) — ai^2(A)e~'*/^ with d the distance between the left and right Majorana 
end modes. The pre-factors aj(A) depend on the local couplings, including the hopping coefficients, pairings, and so 
on, between electrons belonging to the same (for j = 1) or different (for j = 2) sectors of the time-reversal partners. 
For the case with constant and homogeneous local couplings, we have that their ratio ai/a2 is proportional to the 
ratio of couplings between electrons of the same and different sectors, and is nearly a constant, justifying the adiabatic 
condition. The above derivative is clearly confirmed with numerical results in the realistic systems with the presence 
of random on-site disorder scattering, as shown in Fig. |S3| 
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FIG. S3: Adiabatic condition and fermion parity conservation for each sector of time-reversal partners in the presence of 
disorder scattering. The energy splitting E between |00) and jll) (red curves) and the ratio R (blue curves) versus (a) the bulk 
gap which varies by tuning the chemical potential, and (b) the distance between the Majorana end modes. In the numerical 
simulation the random on-site disorder potential is considered, with the potential amplitude Vdis ~ l.OmeV. Other parameters 
in the nanowire are taken that Ap = l.OmeV, As = O.SmeV, and Eso = O.lmeV. The coupling energy E is tuned from to 
l.OmeV within the time 1.0fj,s. 

It is noteworthy that for fixed parameter A, the physics of the fermion parity conservation for each sector can also be 
easily understood. For Dill class topological superconductor, the helical p-wave pairings occur between two electrons 
belonging to the same sector of the time-reversal partners. While the change by one in the fermion number of each 
sector conserves the total fermion parity of the system, it changes fermion parity for each sector, and thus breaks a 
p-wave Cooper pair in each sector. This process costs finite energy and is thus suppressed by the p-wave pairing gap 
if Ap dominates over Ag and the time-reversal symmetry is not broken. The previous study in this section further 
proves this conservation law when the coupling Hamiltonian between Majorana end modes is adjusted adiabatically. 

The fermion parity conservation for each sector has several important physical implications. First, this conservation 
law shows that an isolated Dill class ID Majorana wire should stay in one of the four fermion parity eigenstates 
germinated by non-local complex fermion operators / - and / ■ , as long as time-reversal symmetry is not broken. In 
particular, one can always prepare a nanowire initially in the ground qubit state |00) or |11) by controlling the couplings 
(e.g. through tuning the bulk gap) between Majorana modes at different ends of the wire, and then manipulate the 
states adiabatically. Furthermore, for a single wire, the two topological qubit states with the same total fermion 
parity, e.g. |00) and |11) are physically different states. In the last section of this Supplementary Material, we shall 
show that these states can be distinguished in the experiment. This is explicitly different from the situation in a chiral 
topological superconductor, where two states with the same total fermion parity are not distinguishable. To simplify 
the notations in the further discussion, we relabel the block diagonal Majorana modes 7 ,7- as 7j, 7j. Accordingly, 
the diagonal complex fermion modes are redefined as fj , fj . 



S-3. BRAIDING STATISTICS 



With the fermion parity conservation for each sector we shall show below that the exchange of Majorana zero modes 
in Dill class superconductor is generically equivalent to two independent processes of braiding Majorana fermions of 
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two different sectors, respectively. The physics can be understood in the foUowing way. First of all, since the bulk 
is gapped, braiding adiabatically the Majorana zero modes, e.g. 71,71 and 72,72, does not affect the bulk states. 
Furthermore, assuming that other Majorana zero modes are located far away from the two exchanged pairs, in the 
braiding we only need to consider the evolution of the Majorana zero modes which are braided. Finally, since the 
fermion parity is conserved for each sector, in the braiding the two complex fermion modes /i and /i are always 
decoupled and their dynamics can be derived independently. In this section we study in detail the braiding matrix 
for the Dili class ID topological superconductor. 

A. Degenerate ground state 

We first construct the generic degenerate ground states for the Dili class ID topological superconductor. Consider 
that the ID Majorana wire has 2M pairs of Majorana zero modes 71, 71; 72, 72; ..., and 72M, 72M, with different pairs 
of Majorana zero modes well separated from each other. With these modes we can define the M pairs of complex 
fermion modes by fj = (72^-1 + i72j)/2 and fj = (72^-1 — *72j)/2, with j = 1, ..., M. It follows that 

r-'fjT^f,, T-'f,T = -f,. (S30) 

For the bulk states, we denote the Bogoliubov-de Gennes quasiparticle operators with positive energies by d^ and 
dri, with d,, = T^^drfT- Let |vac) be the vacuum state with respect to the electron operators and we construct the 
wave function 

Here vW is the normalization factor. It is clear that \ip) is an eigenstate of ufj and / /,- with eigenvalue Uj and 
fij, respectively. Note that Uj and fij can be either 1 or 0, and their magnitudes depend on the convention used in 
defining the Majorana wave functions Oil]. To be concrete, we shall use the convention that 

|i^) = |ii...i>|ii...i). (S32) 

The ground state with Uj — fij — is then constructed by 

|00...0)|66...0) = /i.../m/i.../i|V'>, (S33) 

and generically 

M M 



B. General results for the braiding process 

Now we study the general results of the braiding process. In particular, we shall show that braiding two pairs 
of Majorana zero modes 7j,7j and 7j+i,7j+i can be reduced to two independent processes of exchanging 7-,-, 7^+1 
and 7j,7j-|_i, respectively. During the braiding process, the Hamiltonian H{X) generically depends on an adiabatic 
parameter A, and at any fixed A-value the time-reversal symmetry is preserved. There are two different contributions 
which may affect the exchange dynamics. One is the Berry phase effect in the degenerate ground subspace, and 
another is that after the exchange, the original ground state subspace may vary and evolve into new forms. This 
study is similar as that in the chiral topological superconductors O |4] , but with the new ingredients of time-reversal 
symmetry and fermion parity conservation for each sector of the time-reversal partners. 

B.l. Berry phase effect 

In the braiding, the Hamiltonian and the many-body ground state vary with the adiabatic parameter A. To 
determine the Berry phase in the braiding, we calculate the Berry's connection by 

An.mW = ih{ni...nM\{ni..SiM\d\\mi...mM)\'fhi...fhM)- (S35) 
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Since we only braid the two pairs of Majorana zero modes 71,71 and 72,72, we can assume that only these four 



Majorana zero modes are A-dependent, while all other Majorana modes are independent of A [4]. Then, from Eq. (S34) 
we know that 

where 

^™irii,mimi(A) = ih{ni; hi\dx\mi; ifii) 

^ ^niSi,mimi(^) +^nini,mimi(-^) + AiiSi,mimi(-^)- (S37) 

For the last line of the above formula we have that 

<k,™.rn,(A) = *fi(^(A)|[//]i-"M/l]i-"M/i]^-"M/i]i-"^9A|V^(A)), 

4!k,„nr^.W = *^(l-"H)(V'(A)|[/J]i-^M/i^]'-"^(aA/i)[/i]^-"M^(A)), (S38) 

4!k,™,*,(A) = ^Kl - mi)(V.(A)|[/t]i-^^[/t]i-«^[/i]i-™^(a,A)|V.(A)). 

From the former section we have shown that the Fermi parity is conserved for each sector of the time-reversal partners, 
and therefore we have An-^ni.mimi = for rii ^ mi or hi 7^ mi. On the other hand, the continuous variation of A 
cannot change the fermion numbers which are discrete values, which implies that 

9a|V'(A))(x|V(A)), (S39) 

and therefore we have 

^l^k,rn,rnM) = ^nkWSn.mJrnrn, , J = 1,2,3. (S40) 

It is easy to check that ^„\i^j(A) is independent of rii and hi. To calculate ^„%j_,„jmj(A) and ^i%j.m^mj(A), we 
need to examine 9a /i and d\fi, which can be generally decomposed as 

aA/i(A) = u{X)fi{X) + u'{\)fl{X) + «(A)/i(A) + «'(A)/7(A) + 5][a(A)d„(A) + a'(A)4(A)] 

a 

+ 5][6(A)J„(A)+6'(A)dl(A)] (S41) 

a 

dxfiiX) = u{X)h(X) + u'(A)/7(A) + «(A)/i(A) + «'(A)//(A) + 5][fi(A)J„(A) + a'(A)dl(A)] 

+ ^[6(A)d„(A)+6'(A)4(A)]. (S42) 

a 

Fermion modes composed of other Majorana zero states are neglected in the above expansion since they are far away 
from /i and /i. Note that the terms corresponding to bulk quasi-particle operators da and d^ cannot contribute to 
the Berry's connection. Therefore for simplicity we also neglect them in the further discussion. From the time-reversal 
transformation of the two complex fermion operators T~^fiT — ./i,T~^/iT = — /i, we have the following restrictions 
in the coefficients 

u{X) = u*iX), li'(A) = K(A)]*, i;(A) = -v*{X), v'{X) = -K(A)]*. (S43) 

The fermion parity conservation for each sector requires that v{X) = v'{X) = 0. Therefore we have 

dxhiX) ^ u{X)fi{X) + u\X)fl{X), 
dxfi{X) = u*{X)fi{X) + u'*{X)fl{X). 



Substituting the above formulas into the second and third lines in Eqs. (S38) we get then 

^i%i,mimi(A) = i^(l - rni)uiX)Sn,„ijA,^„ 

^i%i.mimi(A) = «^(1 - ml)u*{X)6n^r,^^Sn^rh^■ 



(S44) 
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Furthermore, it can be shown that 

u{X) ^ ii^ixMid^hMiX)) = ^(^(A)|aA(7i72)|V'(A)) 

= 0, (S45) 

where we have used the results that 7^ (A) = 1 and the overlapping between 71 and 72 is negligible. With these results 



1(2,3) 



in mind we have A^^ '-' ^ ^ (A) = 0, and conclude that 

= ih{'lp{X)\dx\tp{X))6nimiSnj^rhiSn2m2Sn2m2---^nMmM^nMmM^ (S46) 

which is independent of rij and fij. Since the Berry's connection is diagonal and identical for all qubit states, the 
Berry phase effect does not bring any nontrivial contribution to the braiding process. This effect is similar as the 
situation in the chiral topological superconductors. 

B.2. Ground state variation 

The nontrivial effect for the exchange of two pairs of Majorana zero modes can be resulted from the fact that each 
ground state itself \ni...nM)\ni---nM) can change after the braiding process. Actually, the final state is generically 
related to the initial one via (after braiding 71,71 and 72,72) 

|ni...nM)|"-l---"-Af)final = 6'**^"^ '"^ Vl'^2-- .njl/) |nin2 •• •"-A/)mitial- (S47) 

We can always choose 9{1, 1) = 0. Since the off-diagonal Berry's connection is zero, the above results show that 
braiding Majorana zero modes 71 , 71 and 72 , 72 is equivalent to two independent processes of exchanging 71 and 
72, and 7i and 72, respectively. For the complex fermion operators, we have /i(Afinai) = e*^'^"'^^"-'/i(Ainitiai) and 
/i(Afinai) = e*^^°^^'/i(Ainitiai)- From the time-reversal symmetry we have that 6(1, 0) = —9(0, 1). This leads to 

\ln2...nM)\ln2...nM) final = |ln2---"M)|l"2---%f)initial, 

|ln2...nM)|0n2...nM)final = e'-^°\in2...nM)\0h2...nM)initia.h 

|0n2---«M)|i^2---%f)final = 6"'"° |07l2 ■ .-"m) 1 1^2---%/)initial, 

|0n2...nM)|0n2...nM)final = \0n2...nM)\0n2-nM)inifmh 

where 9q = 9{1, 0) = —9{0, 1). In the next subsection we shall prove that ^o = '"'/S. Then the braiding matrix is given 
by Ui2 = e^'^iT^e^'^'iTa^ which explicitly respects the time-reversal symmetry. 

C. Braiding phases 

Since the exchange dynamics is proved to be equivalent to two independent processes of braiding the Majorana 
zero modes 71,71 and 72,721 respectively, we can construct a toy model of the Dili class Majorana quantum wire to 
study the braiding phases. The simplest case is a two-copy of ID Kitaev model [5] with spin-1/2 fermions respecting 
the time-reversal symmetry. A T-junction is needed and can be formed by a vertical wire (along y axis) which has 
N sites, and a horizontal wire (along x direction) which has 2N + 1 sites. The A^th site of the vertical wire connects 
to the A^ + 1th site of the vertical wire. Moreover, we consider that in the topological region of the T-junction, the 
hopping term equals to the pairing term t — A. We model the Hamiltonian that 

N 2N+1 

^ = -/^ E 4,<rCy,. + \t\ ^ [(e^cl,^ + e-^c,.^)(e^c^i_^-e-^c,+i,^) 

(T;y— 1 a;— 1 

+ (e"^44 + e*c,4)(e"*4+i.; - e^c.+i4)] . (S48) 

where (f> is the hopping and pairing phase for the vertical wire and /^ < is the chemical potential for the horizontal 
wire. The chemical potential for the horizontal wire, the hopping and pairing in the vertical wire are set to be zero. 
In this configuration the vertical wire of the T-junction is initially in the trivial phase and the horizontal wire in the 
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topological phase. Note the braiding dynamics should be independent of </>, we shall consider (j) — ioT simplicity, 
and then 

N 2JV 

cr;y— 1 x—1 

With the above model, we have four decoupled Majorana zero bound states, 71,^,^ and 72Ar+i,t,^; localized at two 
end sites. The initial four degenerate ground states are given by 



11)11) 
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litial 



itial 
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22N 
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0"=t:i p=Qli<...<l2p 
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p=0 !i<...<J2p 

N 2N+1 



'2p,<T ■■ 2l,0" 



^i2p,f'^il:t 



1 



Z^ Z^ ^p+uf-^n.-t -"- + 2^ Z^ "^i; 



p=0ii<...<i2p+l 

N 2N+1 



N 2N+1 

2^ Z^ ^i2p+l,i'"^il 

p=0ii<...<i2p + l 

N 2N+1 



J 



i2p,l"' il,l 



initial 



22N 11 2-^ 2-^ '^i2p+l'-'^il, 



p=0 ii<...<i2p 



vac) 



I vac). 



It is straightforward to verify that for all above states the average number of electrons is A^ = 2N + 1. 

The Majorana end modes can be transported by tuning adiabatically the parameters fi and t {— A) in the T-junction 
[3]. After the braiding the Majorana modes 7i,-f4 and 72Ar+i,t,i exchange their positions, which can be pictorially 
described by reversing the hopping and pairing direction in the Hamiltonian ( S49 ) . This implies that after braiding 



the new Hamiltonian is obtained by taking t — > — i, since the hopping and pairing to the left and right directions in 
Eq. ( S49 ) explicitly have opposite sign. Therefore we get the new Hamiltonian 



N 2N 

(T;y— 1 x— 1 



(S50) 



The corresponding new ground states then read 

N 2N+1 



|1)|1)™ = 

|i)|6)™ = 



^ n 1+E E M)' 



22N 
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1 + E E (-DHp^f-lr E E (-iKp.. 
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N 2N+1 -| p N 2N+1 
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vac) 



N 2N+1 

|o)i6)™ - ^ n E E (-i)^ip,.-cl. 



vac) 



To finish the braiding, we need to adiabatically transform H' back to the initial form H. This can be performed by 
considering the following Hamiltonian 



N 



2N 



H, 



M / , ^y.ty^y-A '' 7 , 



(7-,y=l 



x=l 



-i\7r/2 J _ iA7r/2 



iA7r/2 t 



(e^^ ^^ 4.t - e^^ c,,t)(e-'^-/^4^i_^ + e'^'^/^^.+i.t) + 



+(e^'"/'cl.; " e-^-/^c..;)(e*^-/2cl+,_^ + e-'^^/^c.+i.;) 



(S51) 
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where A is an adiabatic parameter. The adiabatic ground states are given by 

N 2N+1 



11)11) 



|i)|o)a = 



IDa = 



A — 



22N 
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^ n 1+E E i~^re- 
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N 2N+1 



-ia^Xirp J t 
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*2p,T *i^T 
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From the former subsection we know already that the Berry's phase 9b is the same for all states. With the above 
states we can check directly that 0b = for all states. The vanishing Berry's phase is because the time-reversal partners 
in each ground state contributes oppositely to the Berry's phase. This is reasonable, since a nonzero diagonal Berry's 
phase actually breaks time-reversal symmetry. For A = we have H{X) = H' and for A = 1 the Hamiltonian 
transforms back to the initial one H{\ ^ 1) ~ H. Therefore at A = 1 we obtain the final ground states by 



1 1) 1 1) final = 1 1) 1 1) initial, 

|l)|0)final = « 1 1)1 6) initial, 

|0)|i)final = -i|0)|i)initial, 

|0)|0)final = |0)|0)i„itial. 

We therefore complete the proof that the braiding phase 6*0 = 0{l, 0) = —0{Q, 1) = 7r/2. 



(S52) 
(S53) 
(S54) 
(S55) 



D. Braiding matrix and applications 



According to the results in Eqs. (S52) to (S55l, the braiding matrix for exchanging 7j,7j and 7^+1,7^+1 can be 
constructed by Ujj+i{T,T) = e^^^'^^+'^e^^^'^^+^, which explicitly respects the time-reversal symmetry. To visualize 
the non-Abelian statistics, we consider now two Dili Majorana chains, with the two pairs of Majorana zero modes 
71, 71 and 72,72 localized in the first chain, and the other two pairs 73,73 and 74,74 in the second chain. The four 
complex fermion modes are defined by 



/i == 2''^i+n2)j2 = 2(73 + «74),/i = 2(71 -«72),/2 



1 



(73 -174). 



(S56) 



They satisfy the relation T ^/i,27~ = /i,2 and T ^/i,27~ = —fi.2- In terms of the complex fermion modes, the 
transformation matrix for braiding 72,72 and 73,73 takes the form 



ft ft 



C/23(T, T) = -(1 + ifUl - iflh + ^hfl - */2/l)(l - iflfl + iflh - ^f2fl + */2/l) 



tf~t 



(S57) 



The Hilbert space of the four complex fermions is spanned by sixteen topological qubit states \'nin2)\nin2) — 
(/i)"'(/2)"'(/i)"'(/l)"'|00)|66), where nj,nj = 0,1. The bases can be explicitly written down in the form 
(|00),/it|00),/l|00),/iV2^|00)|) (g> (|06),/i^|00),/2^|00),/7/2^|00)|). With this basis we have further 



U23{T,T) = - 



1 -J1 




ri q 


1 -i 




1 i 


-i I 


® 


i 1 


-i 1 . 




_i 1_ 



(S58) 



Using the above braiding matrix and for an arbitrary initial state we can obtain the final state straightforwardly. If 
the initial state is |00)|00), for instance, we get 



f/23(r,f)|oo)|oo) = ^(|oo)|oo) + |ii)|ii) + i 



;(|06)l|06)j 



ii)-i|ii)|o6)) 
ii)L\n)R + i\io)L\io)R-i\oi)L\oi)R), 



(S59) 
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where the indices L and R represents the left and right Majorana chains, respectively. It is interesting that the above 
state is generically a four-particle entangled state, which carries rich quantum information depending on different 
measurement strategies (see the discussion on measurement in the next section). First, for each Majorana wire if we 
do not distinguish the states of the same total parity (e.g. |10)l and |01)l) in the measurement, in the right hand side 
of the above state the former two terms are equivalent, and can be denoted as |icvon)|fiovon), which implies that both 
Majorana wires are in the even parity state. Similarly, the later two terms are also equivalent, and can be denoted by 
l^odd) |^odd)5 implying that both Majorana wires are in the odd parity state. With these notations one can reduce 
the original state into an effective two-qubit entangled one 

t/23(r,f)|00)|00) = -^(|Leve„)|i?cvc„) + l^odd) |i?odd)) , (S60) 



On the other hand, to identify the state (S57) as a four-particle entangled one, we should be able to measure the 
fcrmion parity for each sector of the time-reversal partners in a single wire. A novel scheme for the measurement 
will be proposed and studied in the next section. Similarly, with three Majorana wires of Dili class, we can generate 
a six-qubit code entanglement through two braiding processes. This shows the natural advantage in generating 
multi-particle entangled state using Dili class topological superconductors, which can be very useful in the quantum 
information processing. For example, a five-qubit code entanglement is the minimum requirement to realize an error 
correcting code [7] . 

Furthermore, a full braiding, i.e. braiding twice 72,72 and 73,73 yields the final state by 

f/|3(T,f)|oo)|66) = |ii)|ii), (S6I) 

which distinguishes from the initial state in that each copy of the p-wave superconductor changes fermion parity. In 
contrast, a full braiding of two pairs of Majorana fermions in a chiral topological superconductor transforms the state 
back to the initial one, which is therefore always trivial. From these discussions we find that in the braiding the 
Majorana modes 7j do not feel their time-reversal partners 7^, which is an essential difference from the situation in 
exchanging two pairs of Majorana fermions in a chiral superconductor, and makes the braiding operator in a TRI 
topological superconductor nontrivial. 

S-4. JOSEPHSON EFFECT IN THE Dili CLASS ID TOPOLOGICAL SUPERCONDUCTOR 

Now we study how to measure the topological qubit states with the Josephson effect. It has been predicted that 
in the chiral ID topological superconductor the Josephson current has 47r periodicity 5j, and the topological qubit 
states for a single wire, |0) and |1), can be detected by measuring the direction of the Josephson currents in the 
junction [HIH]. In this section we predict a novel important effect in the Josephson junction formed by Dili class ID 
topological superconductor, which provides a feasible scheme to detect the topological qubit states in a TRI Majorana 
quantum wire. 

A. Effective coupling Hamiltonian 

We consider a Josephson junction formed by two Majorana nanowire ends with a phase difference = (J)r — 4>l, as 
illustrated in Fig. 3(a), and derive the effective coupling Hamiltonian for the Majorana end modes localized at the 
left (L) and right (i?) ends. The electron tunneling process in the junction is described by 

Ht - Tc[^cfl,i + T£[ jY^flj + H.c, (S62) 

where cl,NtCl.n and cr^i^cr,! represents the electron operators for the A^th site at left and 1st site at right wires 
of the junction, respectively, and T is the tunneling coefficient across the junction. The Majorana end modes can be 
generically expanded in terms of electron operators 

^L = ^{-^LaCLa+ul^^C^^j), ^L^^iuL^jCL^j^-uljC^.^), (S63) 

3 J 

IR = ^{ur.jCrj + U*R,jC^R,j), IR = ^{uRjCR^j + S^j4,i)' (S64) 

j 3 

where U[^/rj = w^ /^ if (pL/R = 0. Note that Cj and Cj represent electron operators of a general time-reversal pair 
at jth site, not necessarily corresponding to spin-up and spin-down, since the spin is not a good quantum number 
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when spin-orbit coupling and s-wave order are present (refer to the discussion in Section II). From the above formulas 
we can solve the electron operators in terms of Majorana and nonzero energy Bogoliubov quasiparticle operators. 
Reexpressing the Bogoliubov quasiparticles in terms of electron operators we can interpret cl MtCl n and cm^cm 
by 



N 



N 



CL,N = ulj^jL - Yl "'LjCL.j - Y^ ^L,jcl,j' 

N N 

CL,N = ul^^jL - Y ^LjCLj - Y '^Ljclj, 

N N 

N N 



CR,i = WR,i7fl - X! "'Rj^Rd - X! ^ 



flj'^-Rj' 



(S65) 
(S66) 
(S67) 
(S68) 



with a constant normalization factor neglected. Here aL/R,jTO,L/B.,j and h]^/iij,bi^/jij are expansion coefficients, 
originated from the quasiparticle operators other than the corresponding Majorana mode. Substituting these results 
into the tunneling Hamiltonian Hx yields that 
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'LJ^LJ 



N N 






i?..4., I+H.c. 



TuL,NU*ii^l{(phLlR + TuL.,NU*R,^ilLlR + H.C 
^ N N s 



- TUl^nIl i J2 ^R^J^Rd + J2 '^kAj ) + "^^fl.i ( J2 "i j^ij + J2 ~^l,3Clj ) iR + H. 
ij(o)+^(i) 



(S69) 



where 



ij(0) ^ 'TuL.Ny'*ii^ilLlR + ^UL,NU*ii^ilLlR + ^-C; 



H^^ 



, N N s . N N 

^i = l J = l ^ ^J = l 3 = 1 



L.J^L^j 



N 



N 



N N s 



In the second equation of the formula (S69) we have neglected the higher-order irrelevant terms. The term 7?'^"^ 



represents the direct coupling between Majorana modes at different junction ends, which gives the first term of the 
effective Hamiltonian Hcs in the main text. This can be seen by noticing that 



UL,N = z|ui,w|e^*-/2, UR,i = \uR,i\e'^'^/^ 

UL,N = -«|uL,Ar|e"^-^/^, Ur^i = \ursW'^'^^'^ , 



with which we can recast _ff ("^ into 



i/(0) = «roCOS-(7L7fl -JlJr), To = 2T|uL,ArU_R,i|. 



(S70) 



(S71) 
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On the other hand, for H^^^ we shall calculate up to the second-order perturbation, which is responsible for the 
second term of ifoff in the main text. From if ^^^ we know that Majorana modes at one end (e.g. the left end) also 
couple to the electron modes at another end (the right end). In the second-order perturbation 7^ and 7l {jr and 
7fl,) couple to cr and cr {cl and cl), respectively. When a nonzero s-wave pairing is present in the quantum wires, 
the electrons ci^/r and cl/r form a Cooper pair and condense. This process leads to an effective coupling between 
Majorana zero modes localized at the same end. Therefore, up to the second-order perturbation in the tunneling 
process, we obtain that 



H. 



(1) -'--^2 

cff =0^ UL,NUL,NlLlL 



+ ^^^Ur^iUr^I'JrJr 



22 "flj"«J 
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H.<S72) 



Here J dT{TT ■ ■ ■) represents time ordered integral. Assuming that the superconducting pairings are uniform in the 
Majorana nanowires, we obtain from the above formula that 
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(S73) 



with £i^/r[As,r; Ap^R]k) the bulk excitation spectra in the left (for L) and right (for R) wires of the junction, 
respectively. The coupling coefficients read 



'^L/r{4') = -'i^^^UL/fi^^N/lUL/R,N/l 



N 
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(S74) 



With the relations obtained in Eqs. (S65) to (S68) we have that aR/LjUR/L^j — \aR/LjaR/Lj\, and bR/L_jbR/Lj 
\bR/L,jbR/L,jW^'*''^^^ ■ Together with the results in Eq. (S70l we can simplify Ti^/ji{(f)) to be 

Tl((/)) = -irie'-^-cc. =risin0, 
Tfl,('/') = -iTie-"f' - c.c. = -Fi sin </), 



and the effective coupling Hamiltonian for Majorana fermions at the same end takes the following form 

H^^ = iVi sin (j){-fLlL - IrIr)- 
The coupling constant Fi is calculated by 

|A.r| 



(S75) 



Fi = T \ul,nul, 
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^s,R\ 



,=1 , £UAs,R;Ap,R;k) 



(S76) 



We have assumed the uniformity of the parameters in the left and right wires of the junction that |As_i| — 
\As,r},\Ap^l\ = |Ap,fl|, and therefore \ul,nUl,n\ = \ur,iUr^i\ and £L{As,L;Ap^L\k) = £r{As,r; Ap^R;k). We note 
that this condition is typically satisfied in the realistic systems. It is clear that the Fi-term vanishes when the s-wave 
pairing A^^l/r is absent. 
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To this end, we combine iJ*^^-' and H^g to reach finahy the effective Hamiltonian for a Josephson junction formed 
by Dili class topological superconductors that 



Hcff = «roCos-(7L7fl 



7L7fl) + *ri sin(/)(7i,7i - 7_R7i^). 



(S77) 



Note that if treating </) as a fixed parameter, the Fi-term in the above formula breaks time-reversal symmetry, which 
also reflects that the leading-order contribution to the coupling between Majorana fermions at the same end (ijjjj) 
should come from the second-order perturbation in the tunneling process. Actually, the direct coupling between jj 
and 7j does not experience the phase difference across the junction and should preserve time-reversal symmetry. This 
is because a uniform pairing phase in one end of the junction can be removed by a constant gauge transformation. 
Therefore the coupling between Majorana fermions at the same end can only be induced by electron tunneling and 
the minimum requirement is to consider the second-order tunneling process. Furthermore, the system restores time- 
reversal symmetry at (j) — rm:, which explains why the Fi-term is proportional to sine/), and has 27r periodicity. 



B. Josephson current 



The Hamiltonian (S77) can be block diagonalized by a constant transformation in the Majorana bases that j[ = 
71+72, 72 = 72 + 71, and 7^ 2 = T"f[^2T~^, which sends Hcs to be iJcff = i(Fo cos 0/2 + Fi sin 0)7^72 -i(Fo cos (/)/2- 
Fi sin 0)7^72. The Andreev bound state spectra are obtained by 



En^,,nf, (0) =(FoCOS0/2 + Fi sin(/))(2n/, - l)+(ro cos 0/2 - Fi sin0)(2nj, - 1), 



(S78) 
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FIG. S4: Joseplison effect in Dili class ID topological superconductor with the inclusion of random disorder scattering, (a) 
The sketch of a Josephson junction with phase difference 0. (b) The single particle Andreev bound state spectra versus the 
phase difference 0. (c) The energy spectra of the four qubit states |nini) (ni,ni = 0, 1) according to the results in (b). (d) 
The Josephson currents (in units of 2eAp//l) for different topological qubit states. In the numerical simulation the amplitude of 
the random on-site disorder potential is set as Vdis ~ l.OmeV. Other parameters are taken that Ap = l.OmeV, A^ = 0.25meV, 
Ebo = O.lmeV, the width of the junction d = 0.75^, and in middle trivial region (gray color) of the junction the chemical 
potential is set to be at the band bottom. 
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which are doubly degenerate at = rmr, reflecting the time-reversal symmetry at these points. Here n ,, ?, are complex 
fermion number operators for / and /' modes, respectively. The Josephson current then reads 

J^ ^{^sm^ + ^cos<j>){2nf, - l)+(^ sin ^ - ^ cos^)(2n^-, ~ 1). (S79) 



The equation ( S79 ) shows all the nontrivial results for the Josephson effect in a Dili class ID topological supercon- 
ductor. In particular, for the even parity states (|06) and |li)) we have the Josephson currents that JS™" — ±f Tq sin |, 
which are of 4tt periodicity, while for the odd parity states (|01) and |10)) the Josephson currents J?'^'* = ±f Ti cos^ 
exhibit 27r periodicity. The difference in the periodicity reflects different mechanisms for the Josephson currents J|™" 
and JS'^'^. The currents JT'^'^ for the even parity states are contributed from the Fo-term in the effective coupling 
Hamiltonian, which is due to the direct coupling between Majorana modes at different ends of the junction. Therefore 
the currents J2^°" are a consequence of the single-electron tunneling process and has An periodicity. On the other 
hand, as contributed from the Fi term, the Josephson currents JT^'^ are resulted from the second-order tunneling 
process which corresponds to the Cooper pair tunneling, therefore being of 2tt periodicity. Furthermore, the currents 
J2 are nonzero even for = 0, which reflects the fact that the odd-parity states violate time-reversal symmetry 
even Hcs preserves at (j) — mn. 

The 47r periodicity of the Josephson currents for even parity states can also be understood in the following way. 
For the direct coupling between Majorana fermions at different ends, when the phase difference across the junction 
advances 2tt, the coupling coefhcients change only tt phase and thus change sign. The generality of this argument 
implies that the An periodicity of JS™" is stable against the disorder scattering without breaking time-reversal sym- 
metry. On the other hand, for odd parity states, the two-fold degeneracy at = mn is protected by time-reversal 
symmetry, which shows that the qualitative properties of the Josephson currents JT^'^ are also stable against the TRI 
disorder scattering. The numerical results are shown in Fig. |S4[ 

With the above results we find that in experiment there are several different ways to distinguish JT™ and JS'^'^. 
For instance, one can measure the periodicity of the Josephson currents, or measure the currents at = n/2 where 
J?™" = ±— I^Fq and JS'^'^ = 0. Furthermore, the two qubit states with the same total parity are distinguished by the 
direction of the currents. The qualitative difference in the Josephson measurements provides direct detection of the 
four topological qubit states in experiment. 
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